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Abstract 

Within a self-consistent proper-time Renormalization Group (RG) 
approach we investigate an effective QGD trace anomaly realization 
with dilatons and determine the finite-temperature behavior of the 
gluon condensate. Fixing the effective model at vanishing temperature 
to the glueball mass and the bag constant a possible gluonic phase 
transition is explored in detail. Within the RG framework the full 
non-truncated dilaton potential analysis is compared with a truncated 
potential version. 

PACS: 05.10.Cc, 64.60.A, 14.70.D 


1 Introduction 

Perturbative QCD studies and lattice simulations demonstrate that the fun¬ 
damental quark and gluon degrees of freedom are the relevant ones at high 
temperatures. Since these degrees of freedom are conhned in the low-tem¬ 
perature regime there must be a quark and/or gluon deconhnement phase 
transition at high temperature. 

Difficulties arise in the description of this phase transition due to the 
absence of local parameters which can be linked to conhnement. QCD at low 
energies can be approximated by effective models which have to satisfy the 
observed symmetry properties and anomaly structure of the theory. Besides 
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the important approximate and spontaneously broken chiral symmetry at 
low energy, the classical QCD-Lagrangian exhibits an additional dilatation 
symmetry or scale invariance in the limit of vanishing current quark masses 
which is embedded in a larger conformal group. This symmetry is broken at 
the quantum level by radiative corrections - the so-called trace anomaly. 

In order to mimic the trace anomaly in an effective framework one intro¬ 
duces a scalar = O’*"’'') dilaton held y with scaling dimension one and a 

logarithmic potential with coupling h of the form [2, ^ 



( 1 ) 


V 


This potential breaks scale invariance and leads therefore to a hnite 
vacuum expectation value (VEV) xo = (0|x|0)- In this way the anomalous 
breaking of scale invariance is effectively taken into account. This can be seen 
by identifying the trace of the energy-momentum tensor, of the effec¬ 

tive theory with the trace of the modihed (symmetric) energy-momentum 
tensor of massless QCD. The trace is related to the divergence of the 

dilatation current which vanishes at tree level where huctuations are 
omitted: 



. ( 2 ) 


Here denotes the non-Abelian held strength tensor, the QCD 

beta (Gell-Mann-Low) function and g is the gauge coupling. 

Beyond tree level, however, massless QCD is no longer scale-invariant. 
Due to the renormalization, a new dimensionful cut-oh hidden in the QCD 
beta function must be introduced. Despite this fact, the anomaly does not 
depend on the chosen regularization scheme. The ehective realization of the 
trace anomaly can be achieved by equating the divergence of the dilatation 
current with the dilaton held. This yields a diherential equation for the a 
priori unknown ehective dilaton potential whose solution is Eq. ©• 

This ehective realization allows for an identihcation of the gluon con¬ 
densate with the VEV Xo of the ehective theory with broken 

scale invariance via Eq. (^). Following the suggestion of Campbell et al. 0 
Xo can be regarded as an ehective order parameter for the deconhning phase 
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transition. The anthors find a strong correlation with the chiral phase transi¬ 
tion, argning that a first-order glnonic phase transition might drive the usnal 
second-order chiral phase transition for two flavors, thereby changing it to 
first order. Their conclnsion is that the qnark condensation temperatnre, 
Tq, for the chiral phase transition is always smaller or eqnal to the glnon 
condensation or confinement temperatnre Tg. Even for large Nc and fixed 
nnmber of fiavonrs the qnark transition wonld be driven by the glnon tran¬ 
sition, becanse Tq and Tg coincide in this case with the critical temperatnre 
of a first-order glnon transition. 

In order to investigate these findings we perform a finite-temperatnre 
analysis of the effective theory within the framework of a proper-time renor¬ 
malization gronp (PTRG) approach. With the logarithmic interaction poten¬ 
tial in Eq. (|l|) the Lagrangian becomes non-renormalizable. One possible 
regnlarization of the infinite T = 0 contribntions wonld be the introdnction 
of a cnt-off parameter A which determines the scale np to which the effec¬ 
tive description of the theory is valid, similar to a model parameterization 
of the Nambn-Jona-Lasinio type. Snch a new dimensionfnl cnt-off wonld 
violate the scaling properties of the effective theory. In order to maintain the 
desired scaling behavior, one has to give the cnt-off-parameter a conformal 
weight of nnit one by some mechanism. 

For the pnrely glnonic Lagrangian specified by the potential in Eq. (]^), a 
direct comparison of the free energy or effective potential at the origin, and 
at the local minimnm, needed in a pertnrbative approach, is not possible 
dne to the breakdown of the Taylor expansion of the potential itself. It is 
therefore impossible to determine pertnrbatively the order (first or second) 
of the phase transition signalled by the disappearance of the scalar dilaton 
field Xo s-s the order parameter. 

In order to address the qnestion of the order of the glnon phase tran¬ 
sition, Campbell et al. 0 introdnced an additional term proportional to 
X” {n < 4) to the effective potential in Eq. (J^) with an a priori nnknown 
temperature-dependent coefficient. They chose phenomenologically a mono- 
tonically increasing function of temperature for this coefficient. Neglecting 
effects of order 0{T^) in the finite-temperature contributions to the effective 
potential, a first-order transition was found, as one would expect from lattice 
simulations [ 3 . 

Within the PTRG approach we can circumvent these difficulties by directly 
calculating the full non-truncated effective dilaton potential for any value of 
the field and not just at the minimum, xo- This in principle allows to analyze 
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the order of the phase transition. 

The outline of the paper is as follows: In Sec. 0 we introduce a self- 
consistent RG method combined with a Schwinger (proper-time) regulariza¬ 
tion. We derive flow equations for two different scenarios. The dilaton model 
is discussed without any polynomial potential truncation in Sec. This full 
potential calculation is compared to an analysis in a truncated potential 
approximation, followed by an error estimation in Sec. Finally, in Sec. § 
we present our conclusions. 


2 The Renormalization Group approach 


The RG approach used in the present work is based on a perturbative one- 
loop expression for the effective potential, which is regularized by Schwinger’s 
proper-time representation of the divergent logarithm. The usual Schwinger 
proper-time integral is modihed by a regulator or blocking function in 
the integrand, thus rendering the resulting flow equation infrared (IR) and 
ultraviolet (UV) hnite. In this way the IR flow scale k is introduced ii- 
The flow equation^ that describes the scale dependence of the action Tfc 
reads[| 


k 


dk 




Tr e 



(3) 


where on the rhs a renormalization group improvement is performed by a 
replacement of the bare second derivative of the action w.r.t. the helds with 
the corresponding running expression T), Q]. This replacement improves the 
one-loop flow equation and corresponds to a resummation of higher Feynman 
graphs i.e. daisy and superdaisy diagrams |^, similar to a Schwinger- 
Dyson resummation. In a plane-wave basis the blocking function serves 
as a momentum regulator similar to the regulator in the ’Exact Renormal¬ 
ization Group’ (ERG) 13 and selects only modes which are peaked around 
the scale k |]^. The direct connection to the cut-off Rk in the ERG for¬ 
malism is still missing but the issue whether this proper-time flow converges 
to the full effective action will not be discussed here. We also omit a detailed 

^For an introductory review see e.g. [Q, Chap. 3 and 4.2 . 

^In general the remaining trace on the rhs of this flow equation runs over all inner 
degrees of freedom. For only one scalar degree of freedom we can omit the trace. 
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analysis of the scheme dependence introdnced by the choice of the regnlator 
fnnction ITT], p^ . 

The regnlator has to fnlhl basically two conditions: The modihed action 
Tfc shonld tend to the fnll action in the IR limit, reqniring that fk^oiT 
cx)) —> 1. In this way the cnt-off is removed and all qnantnm flnctnations are 
taken into acconnt. In addition we have to set /^(r = 0) = 1 for arbitrary 
k. This does not regnlarize the UV regime (r = 0) which is not necessary 
anyway if we start the evolntion at a hnite (large) UV scale A. Finally, for 
the derivation of all flow eqnations we will employ the following choice for 
the blocking fnnction fk in d space-time dimensions: 


k 


dfk 

dk 


4 

dT{d/2) 


(^^2)d/2+lg-rfc2 


(4) 


corresponding to the notation in |^. 

Of conrse, it is not possible to solve the fnll RG eqnation (|^) withont 
trnncation. In order to obtain a tractable set of conpled flow eqnations we 
perform a derivative expansion of the effective action which takes the generic 
form 


= J d'^x\Vk{4>) + ]^Zk{(f)){d^(f)f + .. .| 


(5) 


In this work we consider the local potential approximation (LPA) in which 
the rnnning of the wavefnnction renormalization Zk is neglected {Zk is set to 
nnity) and only the effective potential 14 is taken into acconnt. 

The task is to solve the resnlting conpled non-linear flow eqnation for 
the potential 14 by starting in the nltraviolet with a priori nnknown initial 
conditions and integrate towards zero momentnm with respect to the scale k. 
In this way all qnantnm flnctnations are effectively taken into acconnt. The 
initial conditions at the UV scale are determined in snch a way that calcnlated 
physical qnantities in the IR are matched to predetermined valnes. 

The generalization of onr RG method to hnite temperatnre within the 
Matsnbara formalism is straightforward. Details can be fonnd in Refs. P, 
Ot n. In the zero-momentnm components of the loop integration of Eq. 


one has to introdnce bosonic Matsnbara freqnencies ujn = 27rnT. In the LPA 
the resnlting momentnm integrals can be performed analytically, leading to 
fractional powers in the hnite-temperatnre threshold fnnctions and hence to 
the dimensional rednction phenomenon 


14, 151. 
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3 The full dilaton potential 


The aim of the next two sections is to derive and solve numerically flow 
equations for the gluon condensate from an effective dilaton Lagrangian with 
a logarithmic interaction potential. In this section we consider the full dilaton 
interaction potential without any polynomial potential truncation and will 
compare the results with a truncated potential calculation in the next section. 

As motivated in the introduction we start with the following effective 
Lagrangian in Euclidean spaceQ 



The potential Vq is parameterized by two constants. The depth is deter¬ 
mined by the coupling h and the curvature is given by the minimum xo- 
The minimum is asymmetric and small oscillations about Xo can be inter¬ 
preted as a massive scalar gluonium (glueball) state. Difficulties arise around 
X = 0 since the potential is unstable (actually not dehned as mentioned in 
the introduction). Therefore, a particle interpretation at the origin is also 
no longer possible. For Campbell et al. 0 this phenomenon is an indication 
that glueballs are no longer the relevant degrees of freedom at that point. All 
physical hadrons (if there is any phase transition in the effective model) are 
dissolved into their perturbative gluon content and this phase corresponds 
to a deconhned phase. 

At the minimum xo the potential is equal to Vo(xo) = —h/4. This allows 
for an identihcation of the coupling h with the bag constant B which is the 
negative of the energy density of the vacuum eyac via the relation 


h 

4 


B = Eo(0) - Voixo) 


(7) 


Since we are able to calculate the potential at the origin, this relation can 
be used to determine B in the IR. 

The hrst derivative of the potential w.r.t. the scalar held y, denoted in 
the following with a prime, Vq = dVo/dx, vanishes when evaluated at the 
minimum. The second derivative Vq{xo) = ^/i/xq = yields the tree-level 

^We neglect in this work any quark i.e. mesonic degrees of freedom. 


6 




glueball squared mass, In this work the glueball mass is taken as an 
input. We choose values around = 1.5 GeV which is motivated by recent 
lattice results |]^. The experimental situation is not clear up to now but 
glueball candidates in the mass region of 1.5 ... 1.8 GeV are favored |j^. For 
instance, a possible candidate could be the /o(1500) mesonQ. The value of 
the bag constant B which serves as a second input is not well known but 
ranges between = 0.14... 0.24 GeV. Thus, using e.g. B^^^ = 0.24 GeV 
and = 1.5^ GeV^ at the UV scale (tree level) yields = 0.34 GeV and 


Xo = 0.154 GeV. In Ref. [|I^ it is observed that the temperature scale, where 
thermal excitations become important, is determined by the value of the bag 
constant at T = 0. For = 0.14 GeV, the onset of a significant shift of the 
minimum is seen at T 0.25 GeV while for 5^/"^ = 0.24 GeV this happens 
only for temperatures above T ~ 0.4 GeV. In fact, it was found in ||^ that 
the value for the “critical” temperature is dominated by the value of the 
bag constant. A larger bag constant results in a higher critical temperature, 
although the authors cannot determine a real first- or second-order phase 
transition. From the difference between the RG evolved potential in the IR 
at the true minimum xo and the potential at the origin we can calculate 
the bag constant by means of Eq. (|^ and find a similar dependence of the 
“critical” temperature on the bag constant. 

In order to investigate the behavior of the dilatons in a RG picture we 
apply the method presented in Sec. || to the Lagrangian of Eq. (|]). For the 
dilaton potential in the LPA with d Euclidean space-time dimensions and 
the smearing function]^ we obtain the following flow equation: 


dtVk = 


(47r)V2r(d/2) d (1 + Vk 


( 8 ) 


where t = ln(fc/A) and the tilde indicates the rescaled quantity w.r.t. the IR 
squared scale e.g. In order to solve this coupled equation 

we discretize the field y for a general potential V{x) ^ grid. To close 

the system of equations we need to derive an equation for the first deriva¬ 
tive vi 0. This generates higher derivatives of the potential on the rhs of 
Eq. d^) up to 14 "'. These higher derivatives are then numerically determined 

"^Other possible candidates are the /o(1370) and /o(1710) mesons. Recently, the possi¬ 
bility of mixing of these mesons with nearby qq states has also been discussed. 

®See also for further details and definitions Refs. O, 01. 
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by matching conditions, based on a Taylor-expansion of i4 and VI at inter¬ 
mediate grid points. In this way we obtain a highly coupled closed system of 
flow equations which can be solved numerically with a 5th order Cash-Karp 
Runge-Kutta method • The advantage of closing the system in this way is 
that we do not need higher than third-order derivatives. Since the expansion 
of a logarithmic potential at the origin in a power series leads to singular 
coefficients for powers larger than three this signals the breakdown of a per¬ 
turbative treatment of the model, as already mentioned in the introduction. 
Via the application of the non-perturbative RG method we can circumvent 
these singularities. 

The initial conditions for the coupling h and minimum xo are determined 
at the UV scale A by hts to the glueball mass and the bag constant B at 
the end of the evolution in the IR. We have chosen A = 2.0 GeV, xo = 0.137 
GeV and a coupling h = (0.51GeV)^ and start with a logarithmic dilaton 
potential (tree level). These initial values result after the fc-evolution in a 
bag value of = 0.276 GeV and a glueball ball mass of = 1.502 GeV. 




Figure 1: The /c-evolution of the potential 14 (left panel) versus y and 
the fc-evolution of the normalized minimum Xo(^)/Ao(^ = 0) versus k (right 
panel), = 0.51 GeV, A = 2.0 GeV , Xo(A) = 0.137 GeV). 

Due to the radiative corrections (quantum fluctuations) the shape of the 
effective logarithmic potential is altered when evolving towards the IR as 
shown in Fig. |l| (left panel). In the IR, the potential becomes convex for all 
X values. This is consistent with the dehnition of the scale anomaly. 

The structure of the threshold functions in the flow Eqs. (|^) is very similar 
to those for the 0{N) symmetric model P] . Due to a negative curvature of the 








potential for small x which is typical for broken symmetries, poles can occur 
in the threshold functions. For instance, Vl!/k'^ starts with negative values for 
small field amplitudes y during the /c-evolution. That is the reason why there 
is a correlation of the initial values h, xo and the UV scale A and why they 
cannot be chosen independently. To avoid a pole in the threshold function we 
have to choose initial values in such a way that 14 "/A^ is larger than -1 and I 4 " 
becomes positive for all field amplitudes y during the evolution towards the 
IR. In the right panel of Fig. ^ the fc-evolution of the normalized minimum 
is presented. One recognizes that the evolution almost stops around k ~ 0.4 
GeV. 



Figure 2: The evolved effective potential I4=o versus y for various values of 
the UV cutoff A (hV4 = 0.51 GeV, xo(A) = 0.137 GeV). 

In order to estimate the theoretical uncertainty of these results we show 
in Fig. 1^ the UV cutoff dependence of the evolved effective potential I4=o 
for values A = 2.0 — 2.3 GeV at fixed coupling = 0.51 GeV and a 
fixed minimum xo = 0.137 GeV. The evolved minimum xo of fhe potential 
is shifted towards smaller values for higher UV cutoffs. Due to the longer 
evolution with increasing UV cutoff the absolute depth of the potential is also 
increasing. Nevertheless the change for the bag constant is small (cf. Tab. 
because, for its calculation via Eq. (|^), only the potential difference enters. 
In Tab. the cutoff dependence of the glueball mass is also listed. The 
decrease of the glueball mass with increasing UV cutoff is more sensitive. 

On the other hand if we keep the UV cutoff A fixed and vary the coupling 
h or the minimum xo at the initial UV scale we obtain in both cases a rather 
weak depencence of B and m^. The glueball mass and the bag constant 
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A [GeV] 

bT 4 [GeV] 

[GeV] 

2.0 

0.276 

1.502 

2.1 

0.270 

1.409 

2.2 

0.264 

1.356 

2.3 

0.258 

1.337 


Table 1: UV cutoff dependence of the bag constant B and glueball mass m^. See 
text for details. 


decrease with decreasing h. The variation of xo is anyway restricted to a small 
interval (within 7 % of xo) in order to avoid a pole in the threshold fnnction 
already mentioned above. Within this interval the glneball mass decreases 
within 2% while the bag constant increases within 2 % for increasing xo- 


3.1 Finite-temperature evolution 

Within the Matsnbara formalism, the hnite-temperatnre version of the self- 
consistent flow Eq. (H) is given by the following expression [H, O 


dtVkix) = 


2 r(3/2) 




(4^)(<‘-i)/T(d/2) d + 


(9) 


The non-integer powers in the threshold fnnctions are typical for the 
hnite-temperatnre version within this approach and are governed by the 
choice of the smearing fnnctions fk- At hnite temperatnre we again solve, 
analogons to the zero-temperatnre case, the how eqnation numerically 
by discretizing the x held. We use the same zero-temperature initial con¬ 
ditions at the UV scale A = 2.0 GeV. Note that this choice restricts the 
hnite-temperatnre predictions of our results for higher temperatures. The 
predictive power of the hnite-temperatnre extrapolation is basically deter¬ 
mined by the structure of the threshold functions in the how equations and 
the value of the UV cut-off. A detailed discussion can be found in Ref. lITi 


With A = 2.0 GeV we can ignore the temperature dependence of the initial 
values at the UV scale up to temperatures around T ~ 250 MeV. 

In Fig. ^we show the temperature behavior of the potential in the IR. The 
dashed curve corresponds to the initial bare potential at k = A. After the k 
evolution we obtain the solid curve, labeled by T = 0 MeV, corresponding to 
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Figure 3: The potential Vk=o = 0.51 GeV, A = 2.0 GeV , Xo = 0.137 
GeV ) versus the x field for different temperatures (AT = 200 MeV, see text 
for details). 


the zero temperature IR potential. We repeat this procedure in steps of 200 
MeV. Below temperatures of the order of 200 MeV no signihcant changes in 
the potential shape are observed. 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 



Figure 4: The temperature-dependent normalized minimum Xo{T)/Xoi'^ ~ 
0) of the potential (left panel) and the glueball mass (right panel) for two 
different bag constants (see text for details). 

In analogy to the restoration of spontaneously broken chiral symmetry 
at high temperature, we regard the temperature-dependent minimum of the 
dilaton potential as an order parameter for the deconhning phase transition 
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under the assumption of factorization (see |^). Deconfinement would thus 
be signalled by a transition from a hnite VEV Xo 7 ^ 0 to a vanishing one 
Xo = 0. 


We observe a change of the normalized ’order parameter’ Xq{T)/xq{^) 
with temperature as shown in the left panel of Fig. § for two different values of 
the bag constant B but equal glueball mass rriy. = 1.5 GeV in the IR. For both 
bag constants the condensate stays almost constant up to temperatures of the 
order 200 MeV which is beyond the chiral phase transition temperature for 
two quark flavors |Q. For the smaller bag constant = 0.149 GeV, which 
we obtain for = 0.34 GeV and xo = 0.068 GeV at the UV scale A = 2.0 
GeV, the condensate starts to decrease earlier. Thus, we can qualitatively 
verify the dependence of the “critical temperature” on the bag constant. In 
contrast to Ref. the onset of the shift in the condensate is seen earlier at 
T ~ 190 MeV. It then decreases almost linearly in both cases. The negative 
slope for temperatures beyond 300 MeV does not depend on the chosen value 
for the bag constant. At these temperatures we lose, however, predictive 
power due to the omission of the temperature dependence of the initial values 
at the UV scale (see also Ref. 0). In the right panel the temperature- 


dependent glueball mass is displayed for two different bag constants. The 
mass as function of temperature is independent on the bag constant below 
temperatures of 180 MeV. For the larger bag constant = 0.276 GeV 
it then decreases but increases again for very large temperatures above 800 
MeV. For the smaller bag constant = 0.149 GeV the mass grows almost 
linearly with temperature around T ~ 350 MeV. At very high temperatures 
such a behavior is expected since perturbation theory is applicable (cf. 0 ). 


It turns out that this behavior is almost independent of the initial A, at least 
for temperatures below 250 MeV. 

If we vary the cut-off A between 1.5 and 2 GeV, while tuning the coupling 
h and xo fo keep the glueball mass and the bag constant hxed, we do not 
observe any signihcant changes below temperatures of 200 MeV. 

Another more general way to hnd a possible RG hxed point solution from 
a continuous how equation starts with the rescaled hxed point how equa¬ 
tion 1^. For this purpose we introduce a dimensionless potential U{x) = 
(47r)^A;“^I4(x) and held y = {An)k~^x ia Eq. (g) and set dtU = 0 which 
dehnes the hxed point. In the LPA and d = 4 this yields the how equation 


4U(x)-xR'(x) 


1 1 

2(1 + U"(x)) • 


( 10 ) 
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This equation has a continuum of solutions which depends on the initial 
boundary conditions. For the solution we have to specify two boundary 
conditions. One is fixed by the necessary equation f/^O) = 0 at the critical 
point. If we choose some value for f/(0) we can now numerically integrate 
Eq. (|T^ out to positive fields x- We almost always find a singularity at 
some critical field x = Xc where the numerics breaks down. The result is 
depicted in Fig. ^ where we find two exceptions. The first peak at ?7(0) = 0 
corresponds to a singular solution where U'\0) diverges. Only for the second 
peak at f/(0) = 1/8 the potential seems to exist for all fields y. At this point 
the second derivative vanishes and the potential itself is non-singular and 
constant (trivial) for all fields y. 



Figure 5: The singularity Xc as function of f/(0). 


The peaks are generated by the threshold functions in Eq. (W) and there¬ 
fore depend on the choice of the smearing function. The second solution is 
exactly at f/(0) = 1/8 for a smearing function of the first kind m = 1 §. 
Higher smearing functions will influence the overall factor of the flow equa¬ 
tions moving the second peak towards zero like 2/(m -|- 1)! |^. We conclude 
that no further non-trivial non-singular fixed point arises within the LPA. 
A similar analysis for three space-time dimensions is omited here since the 
extrapolation and implementation of the anomaly is far from trivial. 
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4 The truncated dilaton potential 

In the previous section we have solved the non-truncated flow equation (|[) 
for the full dilaton potential in the LPA by discretizing the scalar held y on a 
grid. Alternatively, one can expand the potential in a Taylor series around the 
non-vanishing local minimum xo (or any other held value) and derive flow 
equations for the (in principle infinitely many) expansion coefficients 
This allows for a direct study of the flow of the coupling constants. We 
will investigate this procedure in this section. Due to the termination of the 
potential series at a finite power we expect new truncation effects. At any 
given finite order of truncation there are always higher-order operators which 
do not evolve and could thus influence the accuracy of the solutions. In the 
limit of infinitely many expansion coefficients we should, of course, reproduce 
all results from the previous section. 

The non-invariant Taylor expansion around the scale-dependent minimum 
Xo (broken phase) up to a specific order M is given by 


22 


Vkix) = Co{k)+ J2^^{x-Xo{k)y 

n=2 


n\ 


( 11 ) 


In the LPA the expansion coefficients, c„, defined at the minimum xo; 


correspond to the n-point proper vertices evaluated at zero momenta ||^ 
Substituting the potential expansion on both sides of Eq. (^ we can deduce a 
set of coupled flow equations for the first several couplings c„, n = 0,2,3,.... 

1 


Co = 


Cl = 


C2 = 


cs = 


C 4 = 


C5 = 


2(47r)2 1 -|- C2/fc2 
Go 


C2X0 + Co 
C3X0 + Co 
C4X0 + Co 
C5X0 + Co 
ceXo + Co 


k^ 

2^ 

k^ 


= 0 

' 


( 12 ) 


^ k^ 


G5 


9 


36- 


.GIGa 




^2 

3CT5 (-T4 


XI 

k<‘ 


- 120 ^ + 240^-60 




Gq 

G,Gl 


90 


k^ 
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The corresponding /3n-functions, 


with the definition 


Gn = 


l^n = ~ dtCn = C, 


(1 + C 2 / k^) 

depend on f3n = /9n(co, Ci,..., Cn+ 2 ) for all n. Due 
to the absence of a reflection symmetry {Z 2 symmetry) we realize that all 
odd and even vertices contribute to the /^n-functions. The expression for Gn 
can be interpreted as a generalized n-point vertex multiplied by a rescaled 
bosonic propagator with the squared dilaton self-energy, C 2 , given at van¬ 
ishing external momenta (cf. [^). This interpretation exemplifies the one- 
loop contributions to the /3„-functions in a transparent and obvious scheme. 
For instance, the first few contributions in the brackets can be visualized 
diagrammatically as 



Figure 6: Diagrammatical contributions to the first several /3„-functions. 

One recognizes the one-loop structure. In the diagrams we have always 
omitted the first term on the rhs of the /9„-functions which is generically 










depicted in Fig. |^. 


Xo 



Figure 7: The term CnXo which has (n — 1) legs. 


The coefficient ci dehnes the minimum xo of the potential and always 
vanishes: 

= c, r 0 , (13) 

We have solved the coupled flow equations (^) numerically with a 5th- 
order Cash-Karp Runge-Kutta method for different truncation orders M in 
order to investigate higher coupling effects in detail. As an example, in 
Fig. H, the fc-evolution of the first coefficients —C 3 , C 4 , the minimum xo 
and the dilaton mass = ^Jo 2 are shown for M = 5. One observes that the 
fc-evolution stops for all dimensionful quantities around the scale k 400 
MeV. All coefficients converge to finite IR values. 

In order to estimate the truncation effects we compare the truncated 
potential calculation for different truncation orders M. It turns out that in 
the truncated potential calculation we cannot choose the same initial values 
as for the full one in the previous sections. With the values A = 2 GeV, = 
0.69 GeV and xo = 0-6 GeV we can indeed perform the truncated potential 
calculation up to the order M = 9 where we have arbitrarily stopped our 
truncation order. The results are presented in Fig. Each curve depicts 
the evolved truncated potential in the IR up to a given order M. Due to 
a slightly different evolution of the minimum xo for different orders M all 
potential curves are plotted versus the x-field normalized in the IR. 

We do not observe an improved convergence of the truncated potential to 
the full potential when increasing the order of the expansion. Obviously, the 
truncation order M = 3 is inappropriate. 

In Fig. 1^ the evolved IR glueball mass, is shown for different truncation 
orders M. All masses are obtained by the same initial values thus displaying 
the influence of higher operators on the mass evolution. We again do not see 
a convergence up to the order M = 9. 
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Figure 8: The scale evolution for the coefficients — Cg' , = C 2 , C3, C4 

and the minimum xo- 


The finite-temperature generalization of the expansion procedure is straight¬ 
forward. In Fig. ^ we show the normalized minimum of the truncated poten¬ 
tial at finite temperature for different values of M. 

For all orders the minimum stays constant up to temperatures around 200 
MeV. For higher temperatures we see truncation effects due to the omission of 
higher operators in the corresponding beta-functions. The rapid temperature 
decrease of the minimum for truncation order M = 3 could also be related 
to the smaller zero-temperature glueball mass at this order. Increasing the 
expansion order M the minimum becomes more stable. Even when varying 
the UV cut-off scale from A = 2 GeV towards A = 1.5 GeV while fixing 
the glueball mass we do not see strong deviations (within 0.1 %) in the 
normalized minima around temperatures of 200 MeV. 

We therefore believe that we can trust the expansion of the dilaton poten- 
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Figure 9: A comparison of the zero-temperature truncated potential calcu¬ 
lation for different expansion order M, demonstrating the influence of higher 
operators to the potential. Each curve is labeled by the corresponding expan¬ 
sion order M. 
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Figure 10; The glueball mass as function of the order, M, of the poly¬ 
nomial expansion. 
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Figure 11: The temperature-dependent normalized minimum /Xo{T = 

0) of the potential for different truncation orders M. (A = 2 GeV, = 0.69 
GeV, xo(A) = 0.6 GeV). 


tial up to a hnite order M when we are interested in temperatures below 200 
MeV. In a forthcoming publication we are going to investigate the hnite- 
temperature behavior of a scaled linear sigma model and study the influence 
of mesons on this result. 


5 Conclusions and Outlook 

We have investigated the thermal behavior of the gluon condensate within 
a self-consistent proper-time regularized Renormalization Group approach. 
An effective realization of the QGD trace anomaly by a scalar dilaton held is 
used, which yields a logarithmic potential parameterized by two constants. 
Due to the breakdown of the effective model for vanishing dilaton held, a 
perturbative analysis is not possible. Therefore, the order of the phase tran¬ 
sition cannot be determine perturbatively. Using non-perturbative RG how 
equations these difficulties can be circumvented and the calculation of the 
ehective potential at the origin becomes possible. 

We have considered two procedures to explore a possible phase transi¬ 
tion in the ehective dilaton model. First the full logarithmic potential is 
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considered. The analysis is then compared to a trnncated potential calcnla- 
tion where trnncation errors are observed and discnssed. The fnll potential 
analysis yields the potential for all dilaton helds and not only aronnd the 
minimnm. In contrast, the trnncated potential calcnlation applies only in 
the vicinity of the global minimnm. 

At vanishing temperatnre, the model has been hxed to the glneball mass 
and bag constant. Onr hnite temperatnre predictions for the fnll poten¬ 
tial are in good agreement with the resnlts of Ref. while the trnncated 


potential calcnlation snffers from slow convergence and artihcial trnncation 
effects. Yet, the hnite temperature results of the truncated potential calcu¬ 
lation correspond qualitatively to those of the full potential analysis. One 
advantage of the truncated version is that a diagrammatical interpretation of 
the beta functions in terms of n-point vertices becomes possible, elucidating 
the underlying physics more clearly. The gluon condensate and the glue- 
ball mass remain unaltered up to temperatures of about 200 MeV where our 
prediction reaches its limit of validity due to the omission of temperature- 
dependent initial boundary conditions in the how equations. 

In contrast to we do not observe any phase transition within the con¬ 
sidered temperature range which makes the use of the expectation value of 
the dilaton held as an order parameter for the gluon deconhnement ques¬ 
tionable. On the other hand this result is not astonishing, since the trace 
anomaly is not expected to vanish at high temperature . 


We have neglected the inhuence of mesonic degrees of freedom which 
should anyhow play a minor role on the gluonic phase transition |T^, 1^, 


The consideration of these degrees of freedom in this RG framework is 
deferred to a forthcoming publication 
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